Gravitational Turbulent Instability of Anti-de Sitter Space 
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Bizon and Rostworowski have recently suggested that anti-de Sitter spacetime might be nonlin- 
early unstable to transfering energy to smaller and smaller scales and eventually forming a small 
black hole. We consider pure gravity with a negative cosmological constant and find strong support 
for this idea. While one can start with a single linearized mode and add higher order corrections to 
construct a nonlinear geon, this is not possible starting with a linear combination of two or more 
modes. One is forced to add higher frequency modes with growing amplitude. The implications of 
this turbulent instability for the dual field theory are discussed. 



Introduction: At the linearized level, anti-de Sitter 
(AdS) spacetime looks just as stable as Minkowski space- 
time or de Sitter spacetime. There are a complete set of 
modes which oscillate and do not grow in time. How- 
ever, at the nonlinear level, it has been suggested that 
AdS might behave very differently. Whereas Minkowski 
space and de Sitter space have been shown to be stable 
under small but finite perturbations [1, 2], this has never 
been shown for AdS. In fact the following intuitive argu- 
ment suggests that it will not be the case. AdS boundary 
conditions act like a confining box. Any finite excita- 
tion which is added to this box might be expected to 
eventually explore all configurations consistent with the 
conserved quantities. This includes small black holes. In 
other words, one might conjecture that after a sufficiently 
long time, any finite excitation of AdS eventually finds 
itself inside its Schwarzschild radius and collapses to a 
black hole [3]. 

In a very interesting paper [4] , Bizon and Rostworowski 
have recently found evidence in favor of this conjecture. 
They considered the spherically symmetric collapse of a 
massless scalar field in AdS. No matter how small they 
made the initial amplitude for the scalar field, their nu- 
merical evolution always produced a black hole. They 
also did a perturbative analysis of the problem and found 
a possible loophole in the above argument. If they started 
with a single linearized mode of the scalar field in AdS, 
they could add nonlinear corrections systematically in a 
way which suggested that no singularity will form. How- 
ever, if they started with a general superposition of lin- 
earized modes and tried to add higher order corrections, 
they had to add higher frequency modes. They also 
showed that under evolution, the energy in the initial 
modes decrease while the energy in the higher frequency 
modes grow. They argued that this is analogous to a 
turbulent instability in which energy is transferred from 
larger to smaller scales. This suggests that a generic fi- 
nite perturbation of AdS will indeed eventually collapse 
to a black hole. 

One drawback of their analysis was that it was re- 
stricted to spherical symmetry, so gravitational degrees 
of freedom were not excited. In addition, one might won- 



der if an angular momentum barrier will prevent a sim- 
ilar collapse at late times. We study the purely gravi- 
tational problem of the nonlinear stability of AdS and 
find results very similar to [4]. Starting with a single 
linearized mode, we compute higher order perturbative 
corrections and find no obstruction. However as soon as 
one has a superposition of linearized gravitational modes, 
higher order corrections force one to add higher frequency 
modes with growing amplitudes. Although we do not fol- 
low the evolution to see whether a black hole eventually 
forms, the following argument suggests that it will. One 
of the singularity theorems shows that closed universes 
are generically singular [5] . This theorem does not apply 
to AdS directly (e.g. because spacelike surfaces are not 
compact) but morally speaking, the negative cosmologi- 
cal constant acts like a confining box for fields inside. So 
one expects that generic solutions will be singular. 

There is no contradiction between this nonlinear in- 
stability of AdS and the positive energy theorem. The 
latter says that if you start with AdS, it will not decay or 
evolve to anything else since it is the unique solution with 
zero energy. The question of whether small amounts of 
energy can collapse to small black holes is very different 
and usually ruled out by arguing that waves disperse at 
late time. This does not happen in AdS. 

The nonlinear generalizations of individual perturba- 
tive modes are geons, i.e., lumps of gravitational energy 
which are held together by their own self gravity. They 
are nonsingular, asymptotically globally AdS, and can be 
viewed as gravitational analogs of boson stars. Since we 
can start with any linearized mode, there are an infinite 
number of geons. However, these solutions are all special 
in that they each have one (helical) Killing field. They 
are exactly periodic in time (answering a question raised 
in [6]). The geons are stable to linearized perturbations, 
but will be unstable at higher order to the turbulent in- 
stability. We will argue that one can put small rotating 
black holes inside the geons and obtain black holes with 
only a single Killing field. (These are purely gravitational 
analogs of the black holes constructed in [7].) Thus, Kerr 
AdS is not the only stationary, asymptotically AdS black 
hole. This had previously been argued for based on the 
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superradiant instability of Kerr AdS. Our small black 
holes sitting inside large geons provide another class of 
examples. 

Physically, one can view the turbulent instability as 
a result of one geon focussing the energy of the other. 
It is perhaps analogous to the collision of plane gravita- 
tional waves which results in singularities in finite time 
due to similar focussing. Of course one difference is that 
plane waves produce singularities after one collision. In 
our case, it takes many collisions of the geons before one 
expects singularities to form. 

Perturbation theory: Wc shall focus on the four- 
dimensional case [8], described by the following action: 




where L is the AdS length scale. We perturb the equa- 
tions of motion derived from (1) by expanding the metric 
about the AdS background, i.e. g — g+J2i h^e 1 , where e 
is a perturbation parameter whose physical meaning will 
be discussed later, and g is the metric of AdS written in 
global coordinates 

g = -(l+£) dt 2 + ^^+r 2 {d6 2 +^ 2 6d<p 2 ). (2) 

At each order in perturbation theory, the Einstein equa- 
tions yield 

Al^=T^, (3) 

where TM is a function of {/i^-* -1 )} and their deriva- 
tives and Al is a second order operator constructed from 
9, 

2Al^S = -V 2 ^-2i?A^S-V a V fc / i «+2V (a V c / i g ) c . 

(4) 

Here, foW = g ab h^, and R a bcd is the AdS Riemann 
tensor. At second order reduces to the familiar 

Landau-Lifshitz pseudotensor [9]. As a consequence of 
the Bianchi identities, VT^j' = for each i. 

We will construct regular finite energy and angular mo- 
mentum solutions of (3) for i > 1. In [10], the SO (3) 
symmetry of AdS was used to show that any regular two- 
tensor, say A, can be expressed as an infinite sum of two 
simple building blocks 

A = E + E + cos ^ sin 0, (5) 

where (A^) represent scalar (vector) type modes, 
which are symmetric two-tensors built from scalar (vec- 
tor) harmonics on the two-sphere. Here (£, m) are the 
usual labels for spherical harmonics, and vector harmon- 
ics on S 2 are of the form *VY^ m . The last term in (5) 
accounts for the fact that we are using a real represen- 
tation for the spherical harmonics, so both the sin <f> and 



cos^ terms need to be included. Each A £ s m (A^ ) 
is parametrized by seven (three) arbitrary functions of t 
and r. 

It is possible to show, by applying the expansion (5) 
to both foW and TW in (3), that any solution of (3) is 
described by two decoupled PDEs of the form 

^l(t,r) + V e W (r)^l(t,r) = T#, (t,r), (6) 

where $e m {t,r) is a Kodama-Ishibashi like (gauge in- 
variant) variable from which h^ m can be recovered (in 
a particular gauge) through a linear differential map 
[10]. One of the equations governs vector-type modes, 
while the other is related to scalar-type modes. Reg- 
ular solutions of (6) are in one to one correspondence 
with smooth solutions of (3). Here, T^L{t,r) is a scalar 
source term, that can be expressed as a function of 
the components of and its derivatives, V^ l \r) is 

a potential that can be found in [10] and D s is the 
d'Alambertian associated with the auxiliary orbit space 
ds 2 = -(1 + r 2 /L 2 )dt 2 + L 2 dr 2 /(L 2 + r 2 ). 

At each order in perturbation theory, we will impose 
regularity of h^ m , seen as a tensor on a fixed AdS back- 
ground. This in turn induces very strong regularity con- 
ditions on any solution of (6). In order for h^ m to have 

a regular centre, ^fL, ~ 0(r e ), as r — >■ 0. At asymptotic 
infinity, we require the metric to be asymptotically glob- 
ally AdS. The asymptotic form of is such that any 
solution of (6) behaves as 

^l-R^ + ^y^ + Oir- 2 ), (7) 

where Re, m and Se tTn are arbitrary functions of t. Note 
that the condition of asymptotically globally AdS is im- 
posed on h^ m not on $^ directly. This means that 

(i) 

we must reconstruct h\ m before imposing the boundary 

conditions on It turns out that if we want our ge- 

ometry to approach the Einstein static universe at the 
boundary, we must keep the leading term and choose 
S e , m (t) = 0. 

A couple of comments about Eq. (6) are in order: i) wc 

(i) 

need to compute all the previous h\ m in order to calcu- 
late the source term in (6) at the next order; ii) if the time 
dependence of the source term is an overall multiplicative 
factor of the form cos(ui£. m t), (6) can be solved by assum- 
ing a separable form for $ty m {t, r) — cos(<D^ m i)i?(r). The 
only exception to (ii) arises if and only if ujg. m coincides 
with any of the AdS gravitational normal frequencies de- 
scribed below. In this case, the modes are said to be 
resonant, and the general solution of (6) now grows with 
time: 

$ ^m(*' r ) = cos (tie,m,t) R Ze,me,rn( r ) 

+ t sm(u) Lm t) 

(r). (8) 
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At first order, ff^ = 0, and one can Fourier trans- 
form Eq. (6) in time, with Fourier parameter uj£, reducing 
the problem to the study of a single ODE of the Sturm- 
Liouville type. For simplicity, we focus on scalar-type 
modes. Requiring a regular centre and the solution to be 
asymptotically globally AdS imposes Luj£ = 1 + I + 2p, 
where p E {0, 1, 2 . . .} is the radial overtone. The positiv- 
ity of ujj indicates that AdS is linearly stable. For p = 0, 
(6) admits a simple solution: 

(i) r e+1 

^,m(*> r ) = — rzrj+rae,m cos(w*t), (9) 

(H + L z ) 2 

where a^ m is a normalization constant. Gravitational 
waves correspond to I > 2. 

Following [4] , we have chosen two types of initial data: 
one is just t — 2, m — 2 and the other is a linear combi- 
nation of t = 2, m = 2 and t = 4, m = 4. For the single 
mode initial data, we choose the normalization and the 
phase so that 

(2L 2 h i t 1 t ) /r 2 )\ r=0 = -9sin 2 0cos - 2<fj . (10) 

We then calculate T^ 2 \ and rewrite it as in (5), with 
a total of six independent terms, corresponding to six 
values for the pair (£, m). This means we have to solve 
for a total of six PDEs of the form (6). The solutions 
can be found and rendered regular both at asymptotic 
infinity and at the centre of AdS. At third order we find 
that T( 3 ) can also be expressed as a sum of six terms, 
five of which behave as the second order terms, but one, 
with l s — m s — 2 and Lw 2 ,2 = 3, is resonant. This could 
lead to the secular behavior exhibited in (8). However 
it turns out that we can set L 3 ^, 2 (r) — by promoting 
the frequency of our initial data to be a function of e 2 : 
Llo 2 = 3 - 14703 e 2 /17920. Thus, to third order, the 
solution is regular everywhere, with no growing modes 
in time. It is invariant under a Killing vector which is a 
slightly shifted version of the symmetry of the linearized 
mode (10): K = d t + fd^. We have computed 
and calculated the energy as a function of the angular 
momentum to fifth order in e: 

_ _ 3 J g ( 4901 J g \ _ 3 / 4901 J g \ 

Eg ~ 2L ^" 7566^)' W2 ~L ^ 3780^' 

where we defined e by J g = ^nL 2 e 2 . It can also be 
checked that this solution obeys the first law of thermo- 
dynamics dE g = ^ dJ g to 0{Jg). From the structure of 
the equations, we expect the situation at higher orders to 
be similar. It is easy to show that there are no resonant 
terms at any even order. We expect one resonant mode 
at each odd order, but by adjusting the frequency one can 
eliminate the growing mode. The result is a (nonlinear) 
geon. 

The two mode initial data behaves quite differently 
At the linear level we start with a composite mode 



$^(t,r) + $£](t,r). At second order we find that 
is now expanded as a sum of 17 terms, none of which are 
resonant. For each of these terms one can calculate the 
right hand side of (6) and a regular solution can always be 
found. The situation changes at third order. Here we find 
that T^ 3 ) can be expressed as a sum of 36 harmonics, 32 
of which behave just as the second order perturbations. 
The remaining four correspond to resonant modes. Two 
of the resonant modes, with Lw 2 ,2 = 3, m s = t s = 2 and 
£0)4,4 = 5, m s = £ s = 4, can be made regular at infinity 
and the origin by a suitable change in the frequencies of 
the initial data, just like we did for the single mode case. 
For the resonant mode with Lwo,o — 1; m s — is = the 
solution is already regular with o(r) = in (8). The 
interesting behavior occurs for the resonant mode with 
the highest possible frequency, L Ci)e t Q = 7, m s = i s = 6. 
In this case the solution can only be made regular if 
£7,6,6 ( r ) m (8) is non-zero, leading to a power law growth 
in time of the initial perturbation, and thus to an insta- 
bility. It is significant that although there are four reso- 
nant modes at third order, only the highest frequency one 
leads to a growing mode. This is the expected behavior 
for a turbulent flow, for which the amplitude of higher 
frequency modes should become larger as time passes by. 
We expect that when this term is large enough, it will 
create other resonant modes with higher frequency and 
that this process will end with the formation of a rotating 
black hole. 

Black holes with only one Killing field: We now 

argue that one can place a small Kerr AdS black hole 
in the core of the geon without disturbing it. Since the 
geon has only one Killing field, the resulting black hole 
will have only a single (helical) Killing field. The argu- 
ment is based on earlier studies where explicit examples 
of scalar hairy black holes were constructed. When the 
black holes are small, their leading order thermodynam- 
ics is accurately reproduced by a non-interacting mixture 
of two components. In one study, the components were 
charged black holes and charged solitons [11]. In another, 
they were rotating black holes and rotating boson stars 
[7] . Similarly, we assume that a small black hole does not 
interact with a large geon. Absence of interaction means 
that the charges E, J of the final black hole are sim- 
ply the sum of the charges of its individual constituents: 
E = Ek+E 9 , J = Jk+Jq- The Kerr component controls 
the entropy and the temperature of the final black hole: 
the geon has no entropy and has undefined temperature. 
Since the geon has only one Killing field and we place a 
Kerr black hole with a Killing horizon at its center, the 
geon's Killing field must coincide with the horizon gen- 
erator of the black hole. That is, the angular velocity 
of the horizon must be fin = This thermodynamic 
equilibrium condition also follows from maximizing the 
entropy for a given total energy and angular momentum. 
Combining this condition with the thermodynamics of 
the two components, it follows that at leading order, the 
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geon carries all the angular momentum of the system. 

Another way to construct black holes with only a sin- 
gle Killing field is through the superradiant instability of 
Kerr AdS black holes with Ltt H > 1 [12]. This insta- 
bility causes certain modes to grow outside the horizon, 
extracting energy and angular momentum from the black 
hole. If you perturb a black hole by a single unstable 
mode, then it will grow and eventually settle down to 
a black hole with "gravitational hair" outside the hori- 
zon. For small black holes, the onset of this instability 
(for a mode with azimuthal quantum number m) can be 
determined by simply setting fijy = ^, where oj is the 
frequency of the linearized geon. (It suffices to consider 
the linearized geon since the "hair" is very weak near 
the onset of the instability.) This condition is motivated 
both by the Killing field argument above and the fact 
that it saturates the condition for a superradiant insta- 
bility: ui < miln- One finds 




This simple argument reproduces the known results for 
black holes with scalar hair. In [7] it was shown that in 
a phase diagram of E vs J, black holes with scalar hair 
exist in the region below the onset of superradiance and 
above the rotating boson star curve. Similarly, we expect 
that vacuum black holes with a single Killing field exist 
between (12) and the geon curve (11). 

Discussion: There is a connection between the turbu- 
lent instability and the superradiant instability. Superra- 
diance causes modes to grow. If you perturb a black hole 
by a single unstable mode, then it will eventually set- 
tle down to a hairy black hole. However, if you perturb 
the black hole with two unstable modes, both modes will 
start to grow. Their interaction will trigger the turbu- 
lent instability which will transfer energy to higher and 
higher frequency. 

Our discussion so far has been entirely classical. Quan- 
tum mechanically, we expect that the transfer of energy 
from large to small scales will still occur, but now will be 
cut off at a frequency equal to the initial energy. So if 
one starts with an initial state consisting of a large num- 
ber of low energy quanta with total energy less than the 
Planck energy, one will never form a black hole. 

We now comment on the dual field theory interpreta- 
tion of our results using gauge/gravity duality. Since we 
have only considered gravity in the bulk, any field theory 
with a gravity dual must exhibit the same turbulent in- 
stability, and transfer energy from large to small scales. 
The fact that this is so universal seems surprising, al- 
though the final outcome of a small black hole can be 
viewed as thermalization in a microcanonical ensemble. 
It is even more surprising considering the fact that our 
dual system is 2 + 1 dimensional, and in this case classi- 
cal turbulence causes energy to flow from small to large 



scales [13]. Our results indicate that (at least at large 
N) strongly coupled 2 + 1 dimensional quantum theories 
in finite volume behave very differently. We have found 
similar results in AdSs, which might appear to contradict 
the idea that the large N limit of super Yang-Mills is in- 
tegrable. However the discussion of integrability does not 
usually include states with energy of order N 2 needed for 
finite backreaction. 

Our nonlinear geons also have a dual interpretation. 
Since the field theory is on a sphere, at zero temperature 
it is in a confined phase in the sense that the free energy 
is order one and not a power of N. A linearized graviton 
can be viewed as a spin two excitation (e.g. a glue-ball). 
The nonlinear geon is dual to a bose condensate of these 
excitations. It is interesting that (at large N and strong 
coupling) these high energy states do not thermalize. 

We have argued that asymptotically AdS solutions are 
generically singular. It would be of great interest to prove 
a singularity theorem establishing this rigorously. 
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